Joint Analysis of Gravitational Lensing, Clustering and Abundance: 
Toward the Unification of Large-Scale Structure Analysis 



Jaiyul Yoo 1 '^ and Uros Seljak 1 ' 2 - 3 ' 4 
1 Institute for Theoretical Physics, University of Zurich, CH -8057 Zurich, Switzerland 
2 Lawrence Berkeley National Laboratory, University of California, Berkeley, CA 94720, USA 
3 Physics Department and Astronomy Department, University of California, Berkeley, CA 94720, USA and 
4 Institute for the Early Universe, Ewha Womans University, 120-750 Seoul, South Korea 

We explore three different methods based on weak lensing to extract cosmological constraints from the large- 
scale structure. In the first approach (method I: Seljak et. al. 2005), small-scale galaxy or cluster lensing mea- 
surements of their halo mass provide a constraint on the halo bias, which can be combined with the large-scale 
galaxy or cluster clustering to measure the dark matter clustering. In the second approach (method II: Baldauf 
et. al. 2010), large-scale galaxy clustering and large-scale galaxy-galaxy lensing each trace the large-scale dark 
matter clustering, and the two can be combined into a direct measurement of the dark matter clustering. These 
two methods can be combined into one method I+II to make use of lensing measurements on all scales. In the 
third approach (method III), we add abundance information to the method I, which is a version of self-calibrated 
cluster abundance method. We explore the statistical power of these three approaches as a function of galaxy 
or cluster luminosity to investigate the optimal mass range for each method and their cosmological constraining 
power. In the case of the SDSS, we find that the three methods give comparable constraints, but not in the same 
mass range: the method II works best for halos of M ~ 1O 13 M0, typical of luminous red galaxies, and the 
methods I and III work best for halos of M ~ 10 14 Mq, typical of low mass clusters. We discuss the robustness 
of each method against various systematics. Furthermore, we extend the analysis to the future large-scale galaxy 
surveys and find that the cluster abundance method is not superior to the combined method I+II, both in terms 
of statistical power and robustness against systematic errors. The cosmic shear-shear correlation analysis in the 
future surveys yields constraints as strong as the combined method, but suffer from additional systematic ef- 
fects. We thus advocate the combined analysis of clustering and lensing (method I+II) as a powerful alternative 
to other large-scale probes. Our analysis provides a guidance to observers planning large-scale galaxy surveys 
such as the DES, Euclid, and the LSST 
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I. INTRODUCTION 

Large-area galaxy surveys like the Sloan Digital Sky Sur- 
vey (SDSS; lllSI]) have enabled high-precision measurements 
of galaxy clustering over a wide range of separation, and 
galaxy clustering has now become a commonly exercised and 
indispensable tool in cosmology. While galaxies are gener- 
ally expected to trace the dark matter distribution up to an 
overall factor on large scales [3], its relation becomes com- 
plicated on small scales, and inferring cosmological parame- 
ters from galaxy clustering measurements are, therefore, ham- 
pered by the complex relation between the galaxy and the 
dark matter distributions, known as galaxy bias b g . In the lin- 
ear regime, measurements of the galaxy power spectrum in 
redshift space [J] provide ways to break the degeneracy by 
constraining the parameter combination / /b g , where / is the 
logarithmic growth rate. However, the recent analysis of the 
SDSS redshift-space distortion yields a relatively large uncer- 
tainties ~ 30% in the linear regime JH], and its full poten- 
tial remains to be realized in future surveys with larger sky 
coverage (see, e.g., JB-O for other approaches to modeling 
redshift-space distortions on small scales). 

Gravitational lensing uses the subtle distortion of back- 
ground source galaxy shapes to statistically map the fore- 
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ground matter distribution that causes the gravitational shear 
(e.g., lfTI14l3tl ). Especially, galaxy-galaxy lensing measures 
the distortion in background source galaxy shapes around 
the foreground lensing galaxies 11441811 . With the statistical 
power present in the SDSS, high precision (20 — 30er) mea- 
surements of galaxy-galaxy lensing signals are typically avail- 
able for various galaxy samples, making it a useful tool for 
cosmology (see, e.g., Ill8l420l0 . Furthermore, spectroscopic 
redshift measurements of the foreground lens galaxies allow 
the galaxy-galaxy lensing signals 7 to be related to the excess 
surface density AE at the lens redshift [21], which in turn is 
related to the galaxy-matter cross-correlation £ ffm . 

The combination of galaxy-galaxy lensing and galaxy clus- 
tering measurements is, therefore, helpful in breaking the de- 
generacy in galaxy bias b g and measuring the matter fluctu- 
ation amplitude as (see, e.g., Il22l429ll ). In this work we per- 
form a systematic investigation of the cosmological constrain- 
ing power that can be derived in the current and future galaxy 
surveys by combining both measurements of galaxy cluster- 
ing and galaxy-galaxy lensing. In particular, we are interested 
in gaining the physical insights of the resulting constraints in 
a model independent way. 

In response to the recent development in numerical simula- 
tions and large-scale galaxy surveys, a halo model based ap- 
proach to modeling galaxy clustering and galaxy-galaxy lens- 
ing has been developed Il22ll32l436ll . The key part of these ap- 
proaches is to assume the halo occupation distribution (HOD) 
and the spatial distribution of satellite galaxies, and to relate 
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TABLE I: SDSS galaxy samples. The approximate numbers for the SDSS galaxy samples are taken l3fl[3lll to represent the SDSS Main, 
LRG, and maxBCG samples. The mean mass M of the galaxy samples is obtained from the galaxy-galaxy lensing measurements I18ll3lll . 
The minimum mass M m i n and the maximum mass M max are obtained by matching the number density and the mean mass of each sample. 
All masses are in units of Ii~ 1 Mq. 
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the dark matter distribution to the galaxy distributions (see, 
e.g., ll37lo . However, many parameters and assumptions of the 
models make it somewhat difficult to untangle the true cosmo- 
logical constraining power. In particular, small-scale galaxy 
clustering information (roughly defined to be below twice the 
virial radius of the largest halos) is used for inferring HOD 
parameters of the model such as the satellite fraction and their 
radial distribution inside halos. While this information pro- 
vides useful constraints on HOD parameters, it is difficult to 
derive any useful cosmological information with it. Moreover, 
the limited number of HOD parameters explored to date may 
artificially provide tighter constraints on the cosmological pa- 
rameters. To avoid this concern we take a simpler approach to 
the problem by approximating galaxy samples as individual 
halos with a certain range of mass. This approach is equiva- 
lent to identifying central galaxies and to removing satellites 
in the galaxy samples. Therefore, there is no useful galaxy 
clustering signal on small scales. 

Without the complication of the spatial galaxy distribu- 
tion and the halo occupation distribution, we can compute the 
galaxy clustering and gravitational lensing signals in a robust 
and model-independent way. Galaxy clustering arises from 
the halo clustering and is tracing a biased version of the dark 
matter clustering on large scales, while galaxy-galaxy lens- 
ing can be split into two regimes. Small-scale galaxy-galaxy 
lensing measures the density profiles of dark matter halos and 
estimates its mass, and large-scale galaxy-galaxy lensing mea- 
sures the cross-correlation of the dark matter and the halo dis- 
tributions. We analyze three different methods to extract cos- 
mological information. In the method I, we use small-scale 
lensing around galaxies to determine the halo mass, which in 
turn determines their large-scale bias using theoretical mass- 
bias relation. Once the bias is known one can use galaxy auto- 
correlation to determine the dark matter clustering. This ap- 
proach was first attempted in [38] using the Main sample of 
the SDSS galaxies with limited success. We will show here 
that higher mass samples offer a better chance of success. 
The method II combines large-scale lensing and clustering to 
eliminate bias, and has been developed in detail in 12711 . The 
method III adds abundance information to the method I, and 



becomes a specific implementation of the cluster abundance 
method lf39l l40h . Traditional cluster abundance methods rank 
clusters by their mass and determine their abundance as a 
function of it. Two main issues in cluster abundance methods 
are getting correct mass of the clusters and determining the 
scatter between the mass observable and the mass, since both 
mass calibration and scatter are completely degenerate with 
cosmological parameters l4lll . One can determine mass cal- 
ibration with small-scale lensing and scatter with large-scale 
clustering analysis, hence the method III uses the same infor- 
mation as the method I with added cluster abundance infor- 
mation. 

In this paper, we investigate the cosmological constraining 
power of combining galaxy clustering and galaxy-galaxy lens- 
ing and focus on finding galaxy samples that are best suited 
for this purpose in the SDSS. Furthermore, we attempt to an- 
swer the same questions in future galaxy surveys by extend- 
ing our analysis to higher redshift. Since numerous ongoing 
and planned future surveys are equipped with deep and wide 
imaging capability, there exists another and potentially more 
powerful way to avoid the complication of galaxy bias and to 
directly map the matter distribution: Using the cosmic shear- 
shear power spectrum ll42l - l44ll . Measurements of the auto- 
correlation of subtle shape distortions, however, have proved 
to be difficult due to the intrinsically weak signal-to-noise ra- 
tio and numerous systematic uncertainties 11451 - 14911 . Since the 
first detections [50-53], only a handful of new measurements 
have appeared, mostly using narrow but deeply imaged areas, 
such as from the Hubble Space Telescope 1 54 15511 . and the 
Canada-France-Hawaii Telescope lf56ll . Repeat imaging of the 
narrow stripe-82 in the SDSS enables measurements of cos- 
mic shear of about 200 square degrees B57I459I1 . and no larger- 
scale weak lensing surveys exist yet. We consider the cosmic 
shear measurements as another component of gravitational 
lensing in future galaxy surveys and compare the resulting 
constraints to the combined constraints of galaxy clustering 
and galaxy-galaxy lensing. In our analysis we do not include 
redshift-space distortions, which is another way to extract the 
information about the galaxy bias. Redshift-space disortions 
suffer from significant nonlinear and scale-dependent bias is- 
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sues i6(ill and need to be understood better before they can be 
used for high precision cosmology. 

The organization of the paper is as follows. In Sec. 
we present our model for the SDSS galaxy samples and dis- 
cuss their physical properties. In Sec. [Til] we present three 
model-independent ways to combine measurements of grav- 
itational lensing and galaxy clustering: Small-scale galaxy- 
galaxy lensing and large-scale galaxy clustering in Sec. IIII Bl 
large-scale galaxy-galaxy lensing and large-scale galaxy clus- 
tering in Sec. MI CI and additional abundance information in 
Sec. HHP I Large-scale galaxy clustering is briefly discussed 
in Sec. IIII Al and the combination of various methods is pre- 
sented in Sec. IIII El In Sec. [IV] we extend our analysis to 
the future galaxy surveys such as the DES, BigBOSS, Eu- 
clid, and the LSST. The cosmic shear constraints are comple- 
mented and compared to the galaxy-galaxy lensing and galaxy 
clustering constraints in the future surveys. We summarize 
our results and discuss the implications for planning large- 
scale galaxy surveys in Sec.[V] Our calculations are performed 
by assuming a flat ACDM universe with the matter density 
Q m = 0.23, the matter fluctuation amplitude as = 0.81, and 
the spectral index n s = 0.97. The matter power spectrum 
shape is kept fixed in cosmological parameter variations. 



II. SDSS GALAXY SAMPLES 

Here we describe our simple model for the SDSS galaxy 
samples, which will be used to investigate their cosmologi- 
cal constraining power. The SDSS |H12l completed its ob- 
servation in 2008 and mapped the sky over 8000 deg 2 in five 
photometric bands to a limiting magnitude r =22.5. We con- 
sider nine galaxy samples that closely match the observed 
SDSS galaxy luminosity-bin samples: The SDSS Main galaxy 
(L1^L6), the luminous red galaxy (LRG), and the maxBCG 
samples (e.g., H [M I61H65IT). However, as the SDSS Main 
L6 sample largely overlaps with the LRG sample, we consider 
only the first five galaxy samples among the SDSS Main sam- 
ples. Furthermore, while there exist faint and bright LRG sub- 
samples [61, 62], we combine both LRG samples into a single 
LRG sample. Finally, the maxBCG samples representing the 
most massive clusters are considered with three subsamples, 
each of which differs in richness threshold and hence in mass. 
Our hypothetical galaxy samples are constructed to represent 
the observed SDSS galaxy samples and cover a wide range of 
mass. The details of the SDSS galaxy samples are described 
inTableU 

We approximate these SDSS luminosity-bin samples as iso- 
lated central galaxies occupying individual dark matter ha- 
los. This approximation is valid for bright galaxy samples, 
as the satellite fraction in typical LRGs is shown to be around 
3—5% lfl9ll66Tl . while the approximation breaks downfor faint 
galaxy samples, where a sizable fraction are satellite galax- 
ies that belong to a group or a cluster of galaxies. However, 
instead of modeling those galaxy samples with more free pa- 
rameters and assumptions, we rely on various methods that 
remove satellite galaxies and identify central galaxies (e.g., 
see |H3l). In this way, we can eliminate the uncertainties as- 



sociated with nonlinear modeling of galaxies on small scales 
and focus on the cosmological constraining power that each 
galaxy sample represents. 

Figure Q] describes our model for the SDSS galaxy sam- 
ples. The top panel shows the relation between the halo 
masses and the central galaxy luminosity. The luminosity 
of the SDSS galaxy samples are obtained by using the K- 
corrected rest-frame ?-band magnitude [68], and their mass 
ranges (Af m j n , Af max ) shown as gray boxes are obtained by 
matching the observed number density n g and the mean mass 
M (thick vertical lines) from the galaxy-galaxy lensing mea- 
surements lfl8l l3lll . The dashed curve shows the best- fit rela- 
tion for the halo mass and the central galaxy luminosity l30ll . 
obtained by analyzing the observed SDSS galaxy clustering 
within the luminosity range shown as the dotted lines. We 
extrapolate the relation (solid) beyond its validity regime to 
represent the LRG and the maxBCG samples, but we modify 
the relation to account for the fact that the central galaxies at 
these clusters are just as bright as the LRGs, while the com- 
bined luminosity of clusters is higher. 

The four bottom panels in Figure Q] show the number den- 
sity n g , the total number AT tot , the average volume V avg , and 
the mean redshift z for the SDSS galaxy samples with their 
mean mass and mass range shown as points and horizontal 
bars, respectively. The number density and the total number 
of galaxies are measured quantities, from which we infer the 
average volume and the mean redshift for each sample. 

Furthermore, we extend our approximation for the galaxy 
samples to all mass range and construct continuous mass-bin 
samples. The continuous mass-bin samples are composed of 
halos with the mass-bin interval A In M = 1.0 at each mass. 
Figure [Hi plots the the number density (solid) of the contin- 
uous mass-bin samples with their mass range (dotted) as a 
function of the mean mass. The mass-bin interval is chosen 
to match the observed SDSS galaxy samples (points). The 
other physical quantities N tot , Vg, vg , and z are shown as the 
solid curves in the other panels. 

Low mass faint galaxies are abundant but probe small vol- 
ume, while the massive luminous galaxies are rare but mea- 
sured at larger distance. Therefore, the total number of galax- 
ies are bounded at very low and high masses. The continuous 
mass-bin samples are good approximations to the observed 
SDSS galaxy samples, with one exception for the LRG sam- 
ple, since the LRG sample is obtained with lower limiting 
magnitude than the Main galaxy samples. We will use the 
continuous mass-bin samples to investigate the cosmological 
constraining power at each mass, but we will compare to the 
nine SDSS mass-bin samples (points) to make connections to 
the observations. 



III. COMBINING GRAVITATIONAL LENSING AND 
GALAXY CLUSTERING 

Large-scale galaxy clustering has been measured with high 
precision, and its theoretical interpretation is simple; it con- 
strains the product of the matter fluctuation amplitude cr 8 and 
the bias of the galaxy sample. In this section, we take this con- 
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FIG. 1: SDSS galaxy and continuous mass-bin samples. We approximate galaxy samples as isolated halos and use the relation between 
central galaxy luminosity and halo mass to match the observed SDSS galaxy samples (discrete points, see Table [I). Top panel: Rest-frame 
r-band magnitude at z = 0.1. Vertical solid lines show the mean mass of the SDSS galaxy samples obtained from the galaxy-galaxy lensing 
measurements The minimum and the maximum masses represented by gray boxes are obtained by matching the mean mass and 

the number density of the SDSS galaxy samples. The dashed curve shows the mass-luminosity relation [30] between the central galaxy 
luminosity and its halo mass, obtained by analyzing the clustering measurements of the SDSS Main galaxy samples. Its validity is limited to 
-22 < M r 01 + 5 log h < -18 (horizontal lines). We modify the mass-luminosity relation shown as the solid curve to implement the LRG 
and the maxBCG samples at high luminosity, but account for the flattening of the luminosity of central galaxies at M > 5 x 10 13 h~ 1 M@. 
Bottom panels: The average volume and the mean redshift of the SDSS galaxy samples are obtained by dividing the observed number density 
by the total number of observed galaxies. For the continuous mass-bin samples (solid), we adopt the mass-bin interval A In M (dotted) that 
closely matches the observed number densities of the SDSS galaxy samples, and we compute the average volume and the mean redshift by 
using the mass-luminosity relation. 



straint from the large-scale clustering measurements of each 
galaxy sample as a base in our cosmological parameter analy- 
sis and combine various gravitational lensing measurements 
to derive further improvements on cosmological parameter 
constraints in a model-independent way. The abundance in- 
formation is also considered in conjunction with its mass mea- 
surements from gravitational lensing. 

To facilitate our understanding of these complementary ap- 
proaches, we split the lensing measurements into two regimes 
(small-scale and large-scale) and call these two different ways 
method I and method II, respectively. We investigate in 
Sees. IIII Bl and IIII CI what information each method can add 
to the cosmological constraining power. The constraints from 



the abundance information (method III) and from the full anal- 
ysis (method IV) of gravitational lensing and galaxy cluster- 
ing are discussed in Sees. IIII Dl and lHI El 



A. Large-Scale Galaxy Clustering 

The signal-to-noise ratio of large-scale galaxy clustering 
measurements can be determined by the number of indepen- 
dent Fourier modes obtainable in a survey volume V s probed 
by each galaxy sample. Using the standard mode counting 
method, we compute the signal-to-noise ratio of large-scale 
clustering measurements, accounting for the shot-noise and 
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the sample variance, as 



dk 4irk 3 V x 
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where /c m ; n = dk = 2tt/V s and a factor two accounts for 
the double counting of Fourier modes due to the reality of 
the galaxy fluctuation field. We modeled the redshift-space 
galaxy power spectrum on large scales using linear theory 
P (fc) = (1 + 2/3/3 + P 2 /5)P g (k) and P g (k) = b 2 g P lin (fc), 
where /3 = f/b g and / is the logarithmic rate of growth. 

Given the mass range (M m i n , M max ) of the SDSS galaxy 
and the continuous mass-bin samples, the galaxy bias factor is 
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and the mean mass of the samples is 

cMm«» 
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where the halo mass function dn / dM and its bias factor b(M) 
are computed by using the Tinker et al. [69] and the Tinker 
et al. [70] fitting formulas, respectively. 

Figure [2] illustrates the galaxy bias and the clustering con- 
straints of the galaxy samples. We use Eq. (O to model the 
galaxy bias factor of each sample. Bias is a monotonic func- 
tion of mass and is rather flat in the low mass regime. While 
the SDSS galaxy samples (points) are at different mean red- 
shifts, their bias factors largely agree with the bias factors of 
the continuous mass-bin samples at z = 0.1, as their redshift 
range is narrow z = 0.1 — 0.3. 

The bottom panel shows the clustering constraint on the 
product C = b g (T R (note Eq. HI is the signal-to-noise ratio on 
square of the product). Fainter galaxies (low mass halos) have 
lower shot-noise due to larger abundances, though the signal C 
is rather fiat due to the constant galaxy bias at low mass. How- 
ever, the volume probed by these faint galaxies is so small 
as seen in Figure Q] that there exist only few Fourier modes 
and hence the clustering constraint in Figure|2]is weak at low 
mass. Especially at low mass, the small volume probed by the 
faint galaxies sets the minimum wavenumber k min close to the 
maximum wavenumber fc max = 0.1/iMpc^ 1 we adopted, and 
thereby the clustering constraints at low mass are discrete and 
weak. 

In the high mass regime, where the galaxy samples corre- 
spond to luminous galaxies and clusters, the galaxy bias factor 
continuously increases with mass, and so does the clustering 
signal C. However, while the volume probed by the clusters 
is assumed constant (Fig. [U, the shot-noise of clusters dra- 
matically increases at high mass due to the exponential nature 
of mass function, degrading the clustering constraints of very 
massive halos. 

Given the volume probed by each galaxy sample and the ra- 
tio of their galaxy power spectrum to the shot-noise, the LRG 
sample yields the best clustering constraint. Dashed and dot- 
ted curves in Figure|2]demonstrates the sensitivity to the max- 
imum wavenumber. For various values of fc max , the trend of 
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FIG. 2: SDSS galaxy bias and clustering constraints. Large-scale 
galaxy clustering measurements constrain the combination of the 
galaxy bias b g and the matter fluctuation amplitude a$. For the fidu- 
cial cosmological model, two panels show the mean halo bias factor 
in Eq. lf2} and the clustering constraints in Eq. {T}, computed at the 
mean redshift for the SDSS galaxy samples (points) in Tableland at 
z = 0.1 for the continuous mass-bin samples (curves). The cluster- 
ing constraints are insensitive to the choice of fc m i n = 2tv/V}^ 3 but 
depend on the adopted value of fc max . Due to small volume probed 
by less luminous galaxies, the clustering constraints are not smooth 
at low mass. 



the clustering constraint in mass remains unchanged. While 
the clustering measurements are more precise in the nonlinear 
regime, its theoretical interpretation becomes more prone to 
systematic errors from nonlinear evolution. 



B. Method I: Small-Scale Galaxy-Galaxy Lensing and 
Large-Scale Clustering 

Galaxy-galaxy lensing measurements on small scales pro- 
vide a robust and a model-independent way to estimate the 
mean mass of the galaxy samples. While the mean matter 
density f2 m may not be well constrained by gravitational lens- 
ing alone, the mean mass is well constrained, as it is propor- 
tional to ui m = il m h 2 . So we model two observable con- 
straints, the large-scale galaxy clustering and the small-scale 
galaxy-galaxy lensing, by using three independent physical 
parameters: two cosmological parameters and the mean mass, 
p = (Q m , cr 8 , M). 

The uncertainty in the mass measurements can be computed 
by adopting the Fisher information technique. Given the fidu- 
cial cosmological parameters, we first compute the total num- 
ber JVtot of foreground lens galaxies shown in Figure Q] that 
will be stacked for galaxy-galaxy lensing measurements. As 
a function of angular separation 8 from the lens galaxies, the 
total number of background source galaxies for shape mea- 
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FIG. 3: Galaxy-galaxy lensing and large-scale clustering constraints, (a) Small-scale lensing measurements constrain the mean mass of the 
continuous mass-bin (solid) and the SDSS galaxy (points) samples. The mass uncertainties are derived by using the lensing measurements 
at projected separation R p = 0.2 ~ 2 h~ 1 Mpc, adjusted for each sample in proportion to (M/Mo) 1 ^ 3 , where Mo is the mean mass of the 
SDSS L4 sample. The dashed line shows the mass uncertainties without the mass-dependent adjustment in R p . (b) Lensing constraints on 
galaxy bias, using the method I (solid & circles) and the method II (dotted & squares). The mass uncertainties are converted into uncertainties 
in galaxy bias (solid) by using the theoretical prediction of galaxy bias b(M). Large-scale lensing measurements (dotted) directly constrain 
the galaxy bias factor of each sample, and the measurements constraints are derived by using the signals at R p = 3 ~ 60ft" 1 Mpc. Remaining 
cosmological parameters are held fixed in both cases, including the mean matter density, (c) Cosmological parameter sensitivity of clustering 
constraint C in each mass, (d) Constraints on the matter fluctuation amplitude from galaxy-galaxy lensing and large-scale galaxy clustering. 



surements is 



and the median redshift can be obtained from 



N bg (9) = JV tot n bg 2ir9A9 



(4) 



where the number density of background source galaxies is 
nbg = 1.2 arcmin~ 2 in the SDSS and we use only those galax- 
ies behind the lens galaxies z s > zi + 0.1, i.e., 

dn bg 



nbg = nbg 



dz 
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The redshift distribution of background source galaxies is nor- 
malized, and we model it with three free parameters (a, 6, c) 
as 



dn, 
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where T(x) is the Gamma function. The mean redshift is 
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by solving the incomplete Gamma function 7(0, x). For the 
galaxy-galaxy lensing measurements in the SDSS, we adopt 
(a, b, c)=(1.34, 2, 1/2) with the mean redshift z — 1.38zo = 
0.42 and the median redshift z m ed — l-31zo = 0.40 from 
& 

The total number of background source galaxies determines 
the uncertainties in shear measurements 
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where the intrinsic shape noise is a- ln t = 0.37 and the respon- 
sivity is 7Z — 1.7 fl8il . Note that these numbers are observ- 
ables and held fixed for cosmological parameter variations. 
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FIG. 4: Indices of the cosmological parameter constraints from three 
different methods of combining gravitational lensing and large-scale 
galaxy clustering (I: small-scale lensing, II: large-scale lensing, III: 
abundance). The clustering constraints in Fig. [2] are included in all 
three methods. Parameter indices are derived by using Eqs. l !17t 
and (24). The value of 7 in the method II depends on the small- 
est scale included in the analysis, and its range is plotted as a gray 
area. 



Therefore, the uncertainties in the mass measurements are 

1 rw) 



M 



dM 



(10) 



and the mean shear signal is computed by using the excess 
surface density AS of the foreground lens galaxy sample as 



1(0) = 



dn bg AE(6>, M) 
where the (comoving) critical surface density is 

„2 1 



4-7rG nris l + zt 



(11) 



(12) 



and (ri, r s , r; s ) are the comoving angular diameter distances 
to the lens, the source, and between the lens and the sources, 
respectively. The factor (1 + zi) arises due to our use of co- 
moving angular diameter distances. The excess surface den- 
sity AS (6, M) on small scales is computed by using the pro- 
jected NFW density profile S(i?) JzElzi as the difference 
AS = £(<i?) — between the mean surface density S 

interior to the disk and the surface density at the same pro- 
jected radius R = 6ri. 

Figure[3]2 shows the uncertainty er ln $ in me mass measure- 
ments of the galaxy samples. We limit the small-scale lens- 
ing measurements to the angle, corresponding to the projected 
separation R p — 0.2 ~ 2/i~ 1 Mpc to avoid complications due 
to the baryonic effects at the center. The angular range is set 
for the SDSS L4 sample and is adjusted for each sample in 



proportion to its virial radius oc Af 1 / 3 . However, the impact 
of this adjustment is small (dashed). The uncertainty in the 
mass measurements is dominated by the total number iV to t of 
lens galaxies and hence the total number Nb g of background 
source galaxies. The LRG sample exhibits a weaker constraint 
cr lnfl 7 than the L5 or the maxBCG samples, slightly deviat- 
ing from the trend shown in the continuous mass-bin samples 
(solid), because the mean redshift of the LRG is higher and 
there exist fewer background source galaxies for shape mea- 
surements. 

Now, in order to combine the mass measurements M with 
the large-scale clustering constraints C and constrain the cos- 
mological parameters, we perform an error propagation anal- 
ysis: 



dC m dC , dC - 
dC = — — dn m + — da a + ^ 77 dM 

Ollm 0(7$ OM 



(13) 



With the matter density fixed (AO m = 0), the constraint on erg 
from the large-scale clustering and the small-scale galaxy- 
galaxy lensing measurements is 



Acre 



dC 
9(Jr 



8M a 



(14) 



and this shows how sensitive the large-scale clustering C is 
to the change in the mean mass of the galaxy sample and the 
matter fluctuation amplitude. 

The sensitivity to halo mass can be read off from Figure [2] 
as dC/dM = <j%(db g /dM) . Since the galaxy bias factor is 
a monotonic function of mass and it is nearly constant at low 
mass end, large uncertainty in mass for faint galaxy samples 
can only contribute little to the uncertainty in the galaxy bias 
factor. Figure [3p shows the uncertainty (solid) in bias predic- 
tion b g constrained by using the small-scale lensing measure- 
ments of the mean mass. An order of magnitude uncertainty 
in mass for the SDSS LI— L3 samples translates into fairly 
good estimates of bias factors, while the galaxy bias factors 
are still better constrained for the luminous galaxy samples. 

However, an accurate estimate of galaxy bias may be irrel- 
evant in constraining the cosmological parameter if the prod- 
uct, the clustering constraint C, is independent of the change in 
the matter fluctuation amplitude. Figure |3p plots the sensitiv- 
ity to the fluctuation amplitude (solid) and to the mean matter 
density (dashed). A high peak of a density field becomes less 
biased as the rms fluctuation amplitude increases. Therefore, 
the sensitivity becomes negative at high mass, and there exists 
a zero-crossing in the derivative dC/d<r$, where the change 
in (Tg is compensated by the change in b g , leaving the cluster- 
ing amplitude unchanged. At this mass the method I cannot 
give any cosmological constraints on the matter fluctuation 
amplitude erg. 

Figure[3]i illustrates the constraint (solid) on erg from com- 
bining the large-scale clustering and the small-scale galaxy- 
galaxy lensing. The insensitivity of the clustering amplitude 
to the matter fluctuation amplitude is reflected as no constraint 
around M ~ 5 X 1O 12 /i _1 M . In [38] this method was ap- 
plied to L4 and L5 galaxies, which we see are close to this 
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zero crossing where no information can be extracted. As a re- 
sult, the derived constraints in [38] are relatively weak. The 
matter fluctuation amplitude can be best constrained by using 
the LRG or the maxBCG samples around M ~ 10 14 /i _1 Af Q . 
At high mass, the low number density n g (and hence the total 
number Ntot) is the dominant source of errors in measure- 
ments of the mean mass and the clustering amplitude of the 
galaxy samples. Since the theoretical prediction of the galaxy 
bias factor is based on lensing mass estimates and the galaxy- 
galaxy lensing measurements are obtained by averaging over 
all galaxies with different large-scale environments, the bias 
prediction is not subject to the sample variance errors and is 
little affected by the systematic errors due to the halo assem- 
bly bias (e.g., ||73rtZ5|])- As long as nearly all halos at a given 
mass are included in the sample, the impact of the halo assem- 
bly bias is small. 

With two observables C and M, we can only constrain two 
parameters. To obtain constraints on the cosmological param- 
eter combination, we define cj = a£ f2^, and the propagation 
equation ( fl"3l is 



where 



dC = 



duj 



dC 



n r 



O m duj dC - 

+ 1^7 dM > 

7 ui dM 



dS} m + — — — da s 

7 <?8 



and the power-law indices satisfy 



2 

a 



dC 



(15) 



(16) 



(17) 



Without loss of generality, we let a = 1, i.e., u> = crgO^. 

Figure |4] shows the power-law index 7 (solid) for the 
method I of combining the large-scale clustering and the 
small-scale lensing measurements. At low mass end, the con- 
straint is insensitive to the mean matter density 7 ~ 0, as 
the clustering amplitude is independent of VL m (shown as the 
dashed line in Fig. |3J:). The index becomes infinity, reflect- 
ing that there is no constraint on the matter fluctuation am- 
plitude around the zero-crossing point. Finally, the combina- 
tion of large-scale galaxy clustering and small-scale lensing 
yields constraints onw = crgf2°; 6 at hi§h mass. The fractional 
uncertainty Aw /w is equivalent to the fractional uncertainty 
Aas/as (Fig.[3]f) if the mean matter density is known, other- 
wise one can think of this method of constraining ui. 



C. Method II: Large-Scale Galaxy-Galaxy Lensing and 
Large-Scale Clustering 

Galaxy-galaxy lensing on large scales provides measure- 
ments of the galaxy-matter cross-correlation, and as shown in 
ll25ll and by combining it with galaxy clustering one constrains 
the parameter combination cj = ogO^ w i m ^ e power-law in- 
dex 7 as a function of scales. On large scales 7 = 1, while if 
one extends the method to nonlinear scales one finds 7 ~ 0.55 



down to 3 — 4/i _1 Mpc ll76Tl . The power-law index of the pa- 
rameter combination cj = ogf^ is thus between 0.55 and 1, 
depending on the smallest scale included in the analysis, as 
shown as the gray area in Figure [4] Here we simply use the 
linear theory for the method II, i.e., 7 = 1, in which no mod- 
eling of the galaxy bias factor is involved (hence the method II 
is not subject to the halo assembly bias). 

Similarly to the small-scale lensing constraint on mass, we 
perform a Fisher matrix calculation to obtain constraints ag, 
where Q — b g £l m a\ is the dependence of the amplitude of 
galaxy-galaxy lensing on cosmological parameters. In this 
case the mean shear signal y(0) in Eq. ( TTTb is obtained by 
computing the excess surface density AS on large scales with 
linear theory cross-correlation £ ffm = b g £ m . The lensing 
signal is summed over the projected separation R p = 3 ~ 
60/i _1 Mpc for all galaxy samples. 

The dotted curve in Figure [3}) shows the constraints on the 
galaxy bias factor b g when the cosmological parameters are 
held fixed (or the constraints on the large-scale galaxy-galaxy 
lensing amplitude Q). The method II works well over a broad 
range of mass, and the constraint from the LRG sample is 
~ 8%, consistent with l76Tl . However, if we limit the lensing 
measurements to R p = 3 ~ 30/i _1 Mpc as the case for the 
SDSS Main samples, the statistical power would be reduced 
by a factor two, and the LRG sample is best suited for the 
method II. 

The large-scale clustering constraints are combined with 
the large-scale galaxy-galaxy lensing constraints by propagat- 
ing errors as 



(18) 



With the matter density fixed (AO TO = 0), the constraint on as 
from the large-scale clustering and the large-scale galaxy- 
galaxy lensing measurements is 
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Figure [3]i presents the constraint (dotted) on erg from combin- 
ing the large-scale clustering and lensing measurements. The 
best constraints on the matter fluctuation amplitude can be ob- 
tained from a fairly broad range in mass, covering the SDSS 
L4 to the maxBCG samples. The method II constraints are 
comparable to the method I, but they differ in the mass scale 
that provides best constraints. 

We note that our analysis includes only statistical errors and 
not systematics (see, e.g., I45l - l49ll ). For example, currently 
the dominating source of error is photo-z calibration and con- 
version from observed ellipticity to shear, which contribute a 
combined 5% error on the amplitude of galaxy-galaxy lens- 
ing signal in the latest SDSS analysis B76I1 . This would sug- 
gest that the amplitude cannot be measured to better than 5% 
accuracy with the current analysis pipeline (already compa- 
rable to the statistical errors in the SDSS), and hence the 
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FIG. 5: Additional abundance constraints. Upper panel shows the 
sample variance in abundance due to the limited volume probed by 
the thresholded samples. Bottom panel shows the constraints on the 
matter fluctuation amplitude by adding additional abundance infor- 
mation. The mean mass and the clustering constraints are combined 
to constrain the minimum mass M m i u , the log-normal scatter <7i n M, 
and the matter fluctuation amplitude as (thick solid). We assume 
finM = 0.5 in the fiducial model. In addition, we consider two sys- 
tematic errors: invisible halos and skewness (S3 7^ 0) in the mass- 
observable relation. If some fraction (a g > 0) of halos are devoid 
of galaxies, they may drop out of the sample, acting as a systematic 
error. Dotted (<J ag = 0.1), dashed (u ag — 0.3), and dot-dashed 
(<rs 3 = 5.0) curves show the degradation in Aa$ due to system- 
atic errors, if we marginalize over those parameters in our modeling. 
The constraint can be greatly improved (thin solid), if the log-normal 
scatter, the invisible fraction, and the skewness are perfectly known 
(A(7i n M = Aq 9 = AS S = 0). 



constraint becomes inflated to Acrg ~ 0.05. The same sys- 
tematic error would also apply to the shear-shear analysis, 
which results in Ao^ ~ 0.05 in the linear regime and some- 
what weaker systematic error Aa$ ~ 0.03 in the quasilinear 
regime, as it scales as <r| lf77ll . However, an additional 5% 
shear-calibration uncertainty in the method I sets a floor at 
~ 7% in mass measurements (since M ~ S 3 ^ 2 ), which de- 
grades the erg -constraint to Acrg ~ 0.04. This suggests that 
method I is less sensitive to the overall calibration error than 
method II and shear-shear method. This may be an important 
consideration for the future surveys too. 



D. Method III: Additional Abundance Information 

Now we add the abundance information in the analysis 
and investigate its cosmological constraining power. Since 
the number density n g of galaxy samples is rather sensi- 
tive to the mass-observable relation, we consider in this case 
thresholded galaxy samples, instead of luminosity-bin sam- 



ples. As our fiducial model, we assume that the thresholded 
galaxy samples are described by two parameters, the mini- 
mum threshold mass M m \ n and the log-normal scatter <7i n M 
of the mass-observable relation, given the cosmological pa- 
rameters: p = (Q m , cr 8 , M min , cti„m)- We set a inM = 0.5 in 
the fiducial model. 

As our observables of the thresholded galaxy samples, we 
consider the abundance n g of the galaxy sample, the mean 
mass M of the galaxy sample from small-scale lensing mea- 
surements, and the large-scale clustering C: O = (n g) M, C). 
The number density of the thresholded galaxy samples can be 
computed, accounting for the scatter in the mass-observable 
relation, as 
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where erfc(a;) is the complementary error function. The mean 
mass and the bias of the thresholded samples are then 
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We use the same mass-observable relation for the central 
galaxies of the thresholded samples as described in Figure Q] 
in computing the volume probed by the thresholded galaxy 
samples. Measurement uncertainties in the mean mass and the 
large-scale clustering amplitude are therefore obtained in the 
same way. While the measurement uncertainty in abundance 
is practically zero in observation, our theoretical prediction n g 
in Eq. (l20l > is based on the globally averaged mass function, 
which differs from the measurements due to the sample vari- 
ance (see, e.g., lT78tl ). Therefore, we treat the sample variance 
in abundance as the measurement uncertainty a rlg in abun- 
dance, and for a galaxy sample covering the volume V the 
sample variance is computed as 



An, 



(R) 



(22) 



where R = (3V A /47r) 1 / 3 and a m (R) is the rms matter fluctua- 
tion smoothed by a top-hat filter with R. 

The upper panel in Figure [5] shows the fluctuation in the 
abundance An g /n g due to the limited volume probed by the 
thresholded galaxy samples. As the volume coverage in- 
creases with the minimum threshold mass (or brighter galax- 
ies), the sample variance decreases accordingly, and the ob- 
served abundance approaches to the global mean value to a 
percent level around M m i„ ~ 10 14 h~ 1 M@. The upturn at 
M > 10 14 /i _1 Af Q reflects the fact that the central galaxies 
of more massive clusters are not brighter, while they are more 
biased. 

Using the matrix form between the observables O and the 
parameters p, a general error propagation analysis can be per- 
formed as 



dO = — dp , 
op 



(23) 
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and assuming AO m = the constraint on erg from the thresh- 
olded galaxy samples can be obtained as 
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The bottom panel in Figure [5] shows the constraint (thick 
solid) on as by adding the abundance information. Due to 
the exponential sensitivity, the derivatives of the abundance 
d \nn g jd In p are larger than the derivatives of the clustering 
amplitude and the mean mass, and the sample variance er- 
ror in the abundance is always an order-of-magnitude smaller 
than the other uncertainties. Consequently, the uncertainties 
in the clustering amplitude and the mean mass dominate the 
error budget in Acg at all mass range, and the uncertainty in 
the matter fluctuation amplitude is set by the uncertainty in 
the clustering amplitude at high mass and the uncertainty in 
the mean mass at low mass. 

Significant improvements (thin solid) can be made in con- 
straining the matter fluctuation amplitude using the abun- 
dance information, if we have extra information about the 
log-normal scatter in the mass-observable relation. At a mass 
of 10 14 /i M© this leads to a factor of two improvement in 
the error. This shows the importance of knowing the scat- 
ter in the cluster abundance method. However, obtaining this 
scatter from other observations is difficult, since other observ- 
ables have different distributions and associated scatters. Our 
approach of extracting it from the clustering is more conser- 
vative, as we discuss in more detail below. Though the con- 
straints are still dominated by the uncertainties in the clus- 
tering amplitude and the mean mass, the sample variance be- 
comes non-negligible in this limit. The power-law index of 
the parameter combination u = <J80 m can be obtained by 
using the matrix form as 



O m D [p = (0^, Mmin, CTlnu)} 

as D [p = (as, M mi „, ai nM )] 



(24) 



and Figure [4] shows the power-law index 7 (dashed) for the 
method III. A unique signature is again there exist a mass 
range, where no constraint on the matter fluctuation ampli- 
tude can be derived (7 ~ 00). Since the change in the matter 



fluctuation amplitude results in not only the increase in the 
overall abundance of halos but also the distortion in the halos 
mass function, there exists a mass range M ~ 10 11 h~ 1 M & , 
at which no constraint on as can be obtained from the abun- 
dance. At high mass, the thresholded samples approach the 
mass-bin samples, and the constraints becomes uj ~ C80° n 6 
as in the method I. 



1. Systematic Errors: Invisible Halos and Skewness in PDF 

Although the abundance information is a powerful tool to 
probe cosmology due to its exponential sensitivity at high 
mass, it is equally sensitive to systematic errors in theoreti- 
cal modeling of the abundance. Here we consider two simple 
cases for potential systematic errors: Invisible halos and skew- 
ness in the probability distribution of the mass-observable re- 
lation. Invisible halos would appear to be a killer for this 
method, since with lensing and abundance information one 
would never be able to tell that some clusters are not present in 
the sample. However, with clustering information one can tell 
the difference. For example, if we assume for simplicity that 
there is no scatter then in a given model at a given abundance 
we measure all of the most massive halos at that abundance. 
If a fraction of these massive halos are dark, they are not in 
the catalog and need to be replaced with less massive halos 
to reach the same abundance. But these less massive halos 
also have lower bias, and the two cases can be distinguished 
if we include clustering information in addition to lensing and 
abundance. In this sense dark halos act in the same way as the 
scatter between mass and observable. Our goal in this sub- 
section is to explore this in more detail to see if there is any 
residual difference between the two. 

The first case is that some fraction a g of halos at each 
mass fail to form galaxies and simply drop out of the observed 
galaxy samples. These invisible halos would act as a system- 
atic error in our cosmological parameter analysis. In the pres- 
ence of the invisible fraction, the halo mass function should 
be modified as (1 — a g )(dn/ dM) in computing the number 
density in Eq. d20"i >. while the mean mass and the bias factors 
in Eq. ( TSTT i remain unaffected. The effect of the invisible frac- 
tion on the mean mass and the bias factor arises solely from 
the change in M mm , if the number density is held fixed. 

The second case of systematic errors is a non-Gaussian 
probability distribution in the mass-observable relation. The 
thresholded galaxy samples are modeled by using the thresh- 
old mass Mmin and the log-scatter a\ n M, assuming that the 
scatter between the observed mass and the true mass is a Gaus- 
sian. Now we consider a deviation from the Gaussian assump- 
tion. A non-Gaussian probability distribution function with 
known cumulants n n can be constructed by using a Gaussian 
distribution as 17911 
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where the mean /1 and the variance n 2 of the Gaussian distri- 
bution are set equal to those in the non-Gaussian distribution. 
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FIG. 6: Degeneracy in additional parameters of the thresholded samples, (a) The fractional change in the mean mass M due to the change 
in the log-normal scatter <7i n M- The minimum threshold mass Af m i n is adjusted to compensate for the resulting change in abundance in all 
panels, i.e., n g is fixed as the value in the fiducial model. (6) The change in the invisible halo fraction a g required to match AM in Panel (a). 
While negative values of a g are plotted for illustration purposes, only non-negative values of a g have physical meaning, (c) The change in the 
skewness S3 required to match AM in Panel (a). The abundance is again fixed as in the fiducial model (S3 = 0). (d) The fractional change in 
the galaxy bias factor due to the change A<ri n M, Aa g , and AS3 in each panel. The invisible fraction parameter a g is largely degenerate with 
the log-normal scatter o"i n M, while the skewness S3 is only degenerate below Afmin < 10 14 /i _1 Mq. 



Using the Edgeworth expansion, we only keep the leading or- 
der correction to the Gaussian distribution, 
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where the third-order Hermite polynomial is H3 (x) = x 3 — 3x 
and the skewness is S3 — k^/k^- In the presence of skewness 
(S3 ^ 0), the change in the number density is therefore 
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and similar calculations can be performed for the mean mass 
and the bias of the thresholded samples. 

For intuitive understanding of the degeneracy between 
these additional parameters and the fiducial model parameters, 
we compute the change in our observables O = (n g , M,C), 



given cosmological parameters and the most well-measured 
abundance fixed. Figure^ shows the fractional change in the 
mean mass M due to the change in log-normal scatter <j\ n M- 
In order to compensate for the change in the abundance due to 
the increase in <7i n M, the minimum threshold mass M m i n has 
to increase. However, the dilution of low mass halos results in 
decrease in the mean mass (thin solid). The impact is smaller 
at low Mmin as the mass function flattens and the mean mass 
is dominated by M S> M m j n . 

Same effects can be achieved by introducing the invisible 
halo fraction a g ^ or changing the skewness of the prob- 
ability distribution S3 7^ 0, without changing the log-normal 
scatter cr\ n M- For a fixed abundance, the invisible halo fraction 

(27) 

a g > (thin solid) in Figure[6f) needs lower minimum thresh- 
old mass and hence provides the same effect of reducing the 
mean mass. Since the magnitude of the change in a g is di- 
rectly proportional to the change in abundance, a fixed A<7i n M 
yields progressively small Aa g at low mass. 

A positive skewness (thin solid) in Figure [6^ puts more 
weight on lower mass and brings more low mass halos into 
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the thresholded sample that mimicks the effect of the log- 
normal scatter on abundance. Similar to the log-normal scat- 
ter, the dilution of low mass halos decreases the mean mass. 
However, since thresholded samples are used, the change in 
the probability distribution affects the samples significantly at 
high mass, but little at low mass. Therefore, the degree of 
change allowed for the skewness can be large |S3<7i n M| > 1, 
especially at low threshold mass. 

Figure [6]i illustrates the level of degeneracy of these addi- 
tional parameters. For a fixed abundance, while the change 
in the mean mass can be masked by the change a g or S3, the 
change in the galaxy bias factor responds in a different way. 
The invisible fraction (dashed) is largely degenerate with the 
change (solid) in the log-normal scatter, such that the invisible 
halo fraction needs to be constrained to the level seen in Fig- 
ure |6t». While the skewness (dotted) in the probability distri- 
bution yields somewhat different effects, the hyper-sensitivity 
to S3 at high mass puts a very stringent requirements on S3. 

The impact of these systematic errors on the erg -constraint 
is presented in Figure [5] The fiducial model assumes no in- 
visible halos (a g — 0) or skewness (S3 = 0). However, with 
incomplete knowledge on the invisible halo fraction a g or the 
skewness S3, we marginalize over each of these parameters, 
and the constraint is inflated: a ag =0.1 (dotted), cr a = 0.3 
(dashed), as 3 — 5.0 (dot-dashed). This demonstrates the level 
of constraints on a g and S3 required to avoid a significant 
degradation on Aug. As shown in Figure [6] the impact of the 
skewness is larger at high mass, and the impact of the invisible 
fraction is larger at low mass, but neither are very important 
at this level of precision, although they may become more im- 
portant for the future data sets. Since the method III is essen- 
tially based on the method I, the same argument applies to the 
method EH, regarding the halo assembly bias: provided that all 
halos of a given mass are included in the sample, the system- 
atic errors due to the halo assembly bias should be minimal. 
Moreover, the assembly bias effect is small at high mass. 



E. Combined Analysis of Galaxy-Galaxy Lensing and 
Large-Scale Clustering 

For easy comparison of various methods for combining 
galaxy-galaxy lensing and galaxy clustering, Figure Q sum- 
marizes the constraints on the matter fluctuation amplitude 
derivable from the SDSS galaxy samples. The constraint (thin 
solid) from the method III appears different from that in Fig- 
ure [5] as it is plotted in terms of the mean mass of the galaxy 
samples, rather than the minimum mass. Since the method III 
makes use of the mass measurements in addition to the abun- 
dance information, it is effectively built upon the method I. 
Therefore, the constraints (thin solid) of the method III en- 
compass the constraints (dashed) of the method I, although lit- 
tle difference arises due to the difference between the thresh- 
olded and the mass-bin samples, especially at low mass. 

The dot-dashed curve shows the combined constraints of 
the large-scale galaxy clustering and the galaxy-galaxy lens- 
ing on all scales (method I and II). The constraints are ob- 
tained by adding the constraints from both methods in in- 
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FIG. 7: Constraints on the matter fluctuation amplitude. Various 
curves show as a function of the mean mass of each sample, the 
constraints on Aa$ from three different methods (I, II, III). The dot- 
dashed curve shows the combination of method I and II, representing 
the constraints by combining the gravitational lensing on all scales 
and the large-scale clustering. The thick solid curve (method IV) 
shows the constraints by adding the abundance information to the 
gravitational lensing and the large-scale clustering constraints. 

verse quadrature. In this way, the large-scale clustering con- 
straints in both methods are doubly counted. However, since 
the uncertainty in the clustering measurements is smaller than 
that in the lensing measurements in method II on large scales 
(and vice versa on small scales) and this trend is reversed in 
method I, the double counting of large-scale clustering con- 
straints affects the final constraints little. 

While the method III (solid) provides the best constraints 
among the three different methods, the combination of 
galaxy-galaxy lensing and large-scale clustering measure- 
ments (method I+U: dot-dashed) provides equally strong con- 
straints, while it is less subject to systematic errors that may 
be present in the method III. Finally, all three methods can be 
combined, representing the full constraints (thick solid) by us- 
ing the galaxy-galaxy lensing and large-scale clustering mea- 
surements. Compared to the method I+II, the abundance can 
bring additional information and tighten the constraints, al- 
beit not much, around the mass range from the SDSS L4 to 
the maxBCG samples. 



IV. FUTURE GALAXY SURVEYS 

Here we extend our analysis performed in Sec. [Ill] to fu- 
ture galaxy surveys like the DES 1 , BigBOSS 2 , Euclid 3 , and 



http://www.darkenergysurvey.org 

2 http://bigboss.lbl.gov 

3 http://sci.esa.int/euclid 
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TABLE II: Specifications of the future galaxy surveys adopted for forecasts. The sky coverage f2 SU rvey is in units of square degrees, and the 
clustering measurements are assumed to be at the mean redshift z c \ . The limiting magnitude is mn m in r-band, and the intrinsic shear noise 
is a-y. The mean number nt g of background source galaxies is per arcminute squared, and the redshift distribution parameters (a, 6, c, z mc d) 
of the background source galaxies are described in Eq. {6^. The last column indicates if the surveys are equipped with spectroscopic redshift 
measurements. 



Survey 


^survey 




TOl im 


n bg 


(7- 


a 


b 


c 


•2m cd 


spec-z 


SDSS 


8000 


0.1 


22.5 


1.2 


0.22 


1.34 


2 


1/2 


0.4 


Y 


DES 


5000 


0.5 


24.0 


12.0 


0.20 


2.0 


1.5 


1 


0.8 


N 


DES+BigBOSS 


1000 


0.5 


24.0 


12.0 


0.20 


2.0 


1.5 


1 


0.8 


Y 


Euclid 


20,000 


0.7 


24.0 


40.0 


0.20 


2.0 


1.5 


1 


1.0 


Y 


LSST 


18,000 


0.9 


27.0 


45.0 


0.20 


2.0 


1.5 


1 


1.3 


N 



the LSST , and we forecast constraints derivable from these 
future surveys by combining galaxy clustering and gravita- 
tional lensing measurements. Compared to our previous in- 
vestigation of the SDSS, the key difference is that these galaxy 
surveys have significantly lower threshold in flux, and they 
can observe fainter galaxies, thereby probing larger volume 
at higher redshift and providing higher statistical constraining 
power. 

In addition to these simple scaling changes, two quantita- 
tive differences arise in the future galaxy surveys. Photometric 
surveys like the DES and the LSST will have a limited capa- 
bility to map the galaxy positions in three dimensional space. 
Provided that there exist a small number of galaxies with 
spectroscopic redshift measurements for photometric redshift 
calibration, galaxy-galaxy lensing analysis can be performed 
without much degradation, compared to the surveys with full 
spectroscopic capacity. However, a full three-dimensional 
analysis of galaxy clustering will be unavailable in these pho- 
tometric surveys. Moreover, we consider a hypothetical deep- 
imaging and spectroscopic survey over 1000 deg 2 , which may 
be available from the overlapping region of the DES and the 
BigBOSS. This combination may allow for both gravitational 
lensing and galaxy clustering analyses. The BOSS and the 
BigBOSS alone are not considered here as they are mainly 
spectroscopic surveys without imaging and hence there are 
too few background source galaxies for gravitational lensing 
analysis in these surveys. 

The other difference in the future surveys and possibly the 
major improvement over the SDSS is the availability of the 
cosmic shear measurements. A large number of background 
source galaxies in the future surveys enable high precision 
measurements of the shear-shear auto-correlation, which have 
been measured only in a very limited region of deep imaging 
in the SDSS ll57l - D9ll . In this section, we account for these 
changes by using the galaxy angular power spectrum in com- 
puting the clustering constraints for photometric surveys and 
by comparing the improved constraints in the future surveys to 
those derivable from the cosmic shear measurements. TableHIl 
summarizes the specifications of the future galaxy surveys we 
adopt in this paper. 



4 http://www.lsst.org 



A. Combining Galaxy-Galaxy Lensing and Large-Scale 
Galaxy Clustering 

While it is difficult to know which galaxy samples with 
what properties will be measured in the future surveys, it mat- 
ters little to our present purposes of using them for gravita- 
tional lensing and galaxy clustering measurements. Here we 
approximate the galaxy samples in the future surveys as the 
mass-bin halo samples and compute the physical quantities 
plotted in Figure Q] The number density n g of the galaxy 
samples is computed by adopting the same mass-bin interval 
A In M = 1.0 as in the SDSS, but the resulting number den- 
sity in the future surveys is lower at high mass, as the mass 
function is computed at higher redshift z c \ described in Ta- 
ble HD 

We use the same SDSS relation between the central galaxy 
luminosity and its halo mass as in Figure Q}z (solid) to com- 
pute the mean luminosity distance that these galaxy samples 
can be measured in the future surveys. The mean luminosity 
distances are larger for the galaxy samples in the future sur- 
veys than those with the same absolute luminosities (or the 
same masses) in the SDSS, because the limiting flux fy lm is 
lower in the future surveys. The average volume probed by 
the galaxy samples is therefore V avg (z) oc /^^/(l + z) 3 , 
where the limiting flux is related to the limiting magnitude 
m Hm = — 2.5 log 10 /ii m + m with additive constant m . We 
compute the mean luminosity distances and the average vol- 
umes for each sample in the future surveys by scaling those 
quantities in the SDSS. We can then readily obtain the total 
number N to t and the mean redshift z of the galaxy samples in 
the future surveys. 5 

For photometric surveys like the DES and the LSST, we 
compute the signal-to-noise ratio of the angular clustering 
measurement, similarly to Eq. ([T]!, but accounting for the lack 



5 We note that the limiting magnitudes of the future surveys in TablelTTlrerJ- 
resent those for the photometric imaging, and those for the spectroscopic 
measurements are shallower. However, since they are scaled with the 
SDSS, we suspect no substantial difference in our projection. More impor- 
tantly, our projection for the future surveys should be taken with caution, 
as many uncertain factors can affect the results presented here. 
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FIG. 8: Improvements of combining the large-scale galaxy clustering and the galaxy-galaxy lensing (methods I, II, III, and I+II) on the a$- 
constraints in future galaxy surveys. Each galaxy survey is indicated by various curves: SDSS (solid), DES (dotted), DES+BigBOSS (dashed), 
Euclid (short dot-dashed), and LSST (long dot-dashed). Compared to the SDSS (solid), the major improvements of the future galaxy surveys 
result from the larger survey volume at higher redshift. The specifications of the future galaxy surveys are summarized in Table ITT1 



of three-dimensional information as 



£ (21 + 1) 
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(28) 



where f s ^ y is the fraction of the sky covered by the survey 
and the noise power spectrum is C™ = 47r/ s i cy /7V tot . Adopt- 
ing the Limber approximation and assuming that the galaxy 
samples are uniformly distribution in a narrow redshift range 
z c i ± Az, we compute the galaxy angular power spectrum as 



a _ 



dz 



H(z) ( dn g 
dz 



P 9 (k) 



r2A 7 



Pa k = 



I 



(29) 

where Ar is the width of the redshift bin Az — 0.05. To be 
consistent with our calculations in spectroscopic surveys, the 
range of the angular multipoles in Eq. (f28b is obtained by us- 
ing the Limber relation k — l/r, given the range of wavenum- 
ber (fcmin, fcmax)- In addition, we adopt k niax = 0.15/iMpc -1 
to take advantage of the fact that the structure probed by the 
future surveys remains in more linear regime than that probed 
by the SDSS. 

Figure [8] describes the improvements of the matter fluctu- 
ation constraints over the SDSS by using the three different 



methods of combining the large-scale galaxy clustering and 
the galaxy-galaxy lensing measurements in various future sur- 
veys (different curves). The key improvements result from the 
larger survey volume and the larger number of background 
source galaxies at higher redshift, although the reduced num- 
ber density increases the shot noise contribution. 

Figures [SJi and [Sf) show the constraints on erg by combin- 
ing the large-scale clustering and the galaxy-galaxy lensing 
measurements on small and large scales, respectively. In addi- 
tion to the improvements in the clustering constraints, a larger 
number density of background source galaxies at higher red- 
shift enables higher precision lensing measurements in the fu- 
ture surveys. Indeed, the contribution of the clustering con- 
straints in the method II are weaker than the contribution of 
the large-scale galaxy-galaxy lensing constraints, whereas in 
the SDSS the trend is opposite. However, in the method I 
the contribution of the small-scale galaxy-galaxy lensing con- 
straints is weaker, as the bias factor is highly constrained due 
to its flat nature at low mass and hence the improvements in 
bias estimates are gradual even with substantial improvements 
in mass measurements of the galaxy samples. 

The situation is similar for the method III, shown in Fig- 
ure^, to the case in the SDSS. As more volume is available, 
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source redshift z angular multipole 1 maximum angular multipole 1 



FIG. 9: Constraints on the matter fluctuation amplitude from the cosmic shear measurements in future surveys. Various curves represent 
different galaxy surveys as in Figure[8]and Table Ull SDSS (solid), DES (dotted), Euclid (short dot-dashed), and LSST (long dot-dashed). Left: 
The redshift distribution of the background source galaxies. Middle: The angular convergence power spectrum (thick) and the noise power 
spectrum (thin). Right: The constraints on the matter fluctuation amplitude as a function of the maximum angular multipole. For comparison 
only, we included the SDSS (solid) for the cosmic shear measurements, although it is hard to measure in the SDSS. Horizontal lines represent 
the non-Gaussian contributions of the lensing trispectrum and the sample variance 1 80-82], both of which dominate the error budget on small 
scales. 



the sample variance in abundance decreases even further in the 
future surveys, making the uncertainties in the number den- 
sity nuisance. Consequently, the uncertainties in the cluster- 
ing amplitude and the mean mass are the dominant factors in 
Aos at all mass range, even though the lensing constraints in 
the mass measurements improve more than the clustering con- 
straints, compared to the SDSS. The best constraints on the 
matter fluctuation amplitude in the future surveys are achieved 
at slightly lower mass than in the SDSS, since much fewer 
galaxies at high mass are available at high redshift than at low 
mass, relatively. 

All of the three methods dramatically improve the con- 
straints on the matter fluctuation amplitude over the SDSS, 
promising sub-percent level measurements. While the 
method III has somewhat better statistical power in the SDSS 
than the method I or II, its statistical power is comparable to 
other methods in the future surveys and is not superior to the 
combined method I+II in Figure [8]i While no systematic er- 
rors are considered in Figure[8] they are relatively well known 
for the methods I and II. The systematic errors we considered 
for the method III in Sec. MI DTI are relatively innocuous, but 
the method III may be subject to other unknown systematic 
errors, which may become dominant when the overall con- 
straints reach the sub-percent level in precision. 



B. Cosmic Shear Power Spectrum from Gravitational Lensing 

In addition to the measurements of galaxy clustering and 
galaxy-galaxy lensing, the future galaxy surveys have a large 
number of background galaxies at high redshift enough to pro- 
vide robust measurements of cosmic shear signals. Here we 
compute constraints on the matter fluctuation amplitude by us- 
ing the cosmic shear measurements in the future surveys and 
compare them to the constraints described in Sec. lIV Al 



The shapes of background source galaxies are subtly dis- 
torted by the foreground matter distribution, and the statistical 
measurements of the background galaxy shapes can be used 
to isolate the distortion from their intrinsic shapes. These cos- 
mic shear signals measured by the distortion in shapes are 
represented by the projected matter fluctuation, or the con- 
vergence k. Using the Limber approximation we compute the 
convergence power spectrum as 

Cf = J dz w 2 {z) p m (k = l -;A , (30) 

where the weight function is 



3H$ n 



W(z) = ^ — 



dz r{z s ) 



(31) 



Here we consider different parameters for the background 
source distribution in Eq. © that are more adequate for higher 
redshift surveys than in the SDSS: a = 2 and b = 1.5. The 
median redshifts of the source galaxies in each survey are 
listed in Table [XT] To derive the constraints from the cosmic 
shear measurements, we compute the Fisher matrix J47i[83ll 
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2 [Cf + C™] 2 d Pa d Pl 3 ' 



(32) 



where the noise power spectrum in this case is C" = jf nt /nbg- 
A few caveats are in order, regarding our forecasts of the 
cosmic shear measurements. Equation d32b accounts for the 
galaxy shape noises as the sole source of errors, while it is 
known that there exist numerous systematic errors associated 
with the instruments ll45l l46ll . the photometric redshift mea- 
surements |49[ HH], and the intrinsic alignments [48]. Our 
optimistic assumption should serve as the best possible con- 
straints derivable by using the cosmic shear measurements. 
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Furthermore, while tomographic measurements of the cos- 
mic shear signals at different redshift bins can improve con- 
straints on the dark energy equation-of-state or its time evo- 
lution (e.g., 12310 . we consider only a single redshift bin for 
cosmic shear measurements in the future surveys, since the in- 
formation contents on the matter fluctuation amplitude should 
be the same as those obtained by using multiple tomographic 
bins. 

The left panel of Figure|9]shows the redshift distributions of 
background source galaxies in the future surveys. While the 
redshift distribution parameters (a, b, c)=(2, 1.5, 1) are differ- 
ent from those in the SDSS, the shape itself (various curves) 
is similar to each other, and only the median redshifts dif- 
fer for various future surveys, as listed in Table [II] The mid- 
dle panel shows the convergence power spectrum C ; K (thick) 
and its noise power spectrum C™ (thin) from the cosmic shear 
measurements. With few background galaxies in the SDSS, 
the convergence power spectrum (solid) is dominated by the 
noise power spectrum, already at I > 20, though there ex- 
ist (2Z + 1) modes to be added per each angular multipole. 
The convergence power spectra in the future surveys are larger 
than in the SDSS, since the background source galaxies are at 
higher redshift. The noise power spectra are smaller, simply 
due to larger number density in the future surveys. 

The right panel shows the as -constraints using the cosmic 
shear measurements as a function of the maximum angular 
multipoles. The erg -constraints improve as more modes are 
included, but they saturate, once the noise power spectrum 
overwhelms the convergence power spectrum. In our most 
optimistic consideration, the cosmic shear measurements will 
provide amplitude constraints from the future surveys that are 
well below 0.1%. It would appear that the cosmic shear mea- 
surements in the future surveys trump the previous three meth- 
ods in constraining the matter fluctuation amplitude, if all the 
systematic errors are under control and full statistical powers 
are utilized as we assumed here. This is mostly a consequence 
of the fact that cosmic shear measurements can extract useful 
information on the dark matter clustering well into the nonlin- 
ear regime at high angular multipole I. 

However, there are caveats to this conclusion. One is that 
baryonic effects also change the predictions at high I, and 
these would need to be understood in detail 18411 . Second 
caveat is that there exist a sample variance and a nonlinear 
contribution of the matter trispectrum to the covariance of the 
lensing power spectrum on small scales [80-82]. High degree 
of correlations between modes implies no further information 
can be extracted from the lensing power spectrum on smaller 
scales. In the right panel of Figure|9] the horizontal lines show 
these contributions to the lensing measurement uncertainties 
in the matter fluctuation amplitude in future surveys. We find 
that most of the experiments saturate the useful information, 
long before they run out of signals, and the constraints from 
the cosmic shear measurements are comparable to those from 
the combined method I+II. 



V. DISCUSSION 

Cosmological methods in large-scale structure are tradi- 
tionally divided into weak lensing, galaxy clustering, and 
cluster abundance. While they are usually presented as sep- 
arate techniques, there exists a considerable overlap between 
three methods. For example, the cluster abundance method 
cannot exist without a proper cluster-mass calibration and it 
is widely accepted that mass calibration based on weak lens- 
ing is required [85]. This combines cluster abundance and 
weak lensing. Second, a nuisance parameter in cluster abun- 
dance method is the scatter between cluster observable, such 
as luminosity (in optical, X-ray or SZ), and the cluster mass. 
This scatter can in principle be determined from external data 
ll65i l86ll . but it is not clear that such an approach is reliable. 
A more conservative approach is to determine scatter inter- 
nally from the clustering amplitude [39, 40, 87]. This method 
thus combines lensing, clustering and counting (abundance). 
It is useful to ask the questions such as what is the mass range 
where we get most of the information, how can we combine 
the information from different mass ranges, what happens if 
we drop individual observations etc. The purpose of this paper 
was to address these questions. 

The Sloan Digital Sky Survey (SDSS) has enabled high 
precision measurements of galaxy clustering and gravitational 
lensing for various galaxy samples that cover a wide range of 
halo masses, so we started our analysis with the SDSS. With 
careful modeling of the SDSS galaxy samples, we have ex- 
amined the cosmological constraining power that is contained 
in each galaxy sample and can be derived by combining both 
measurements of galaxy clustering and gravitational lensing. 
Furthermore, we have extended our analysis to other ongo- 
ing and future galaxy surveys with larger volumes at higher 
redshifts such as the DES, Euclid, and the LSST Our joint 
analysis of gravitational lensing and galaxy clustering pro- 
vides a guidance to observational applications and planning 
future large-scale galaxy surveys. However, in doing the anal- 
ysis we had to make a number of simplifying assumptions, so 
our results should be considered as a useful guidance rather 
than conclusive. 

We have investigated the cosmological constraining power 
in three model-independent ways to gain physical insights of 
its information contents, using constraints on the amplitude of 
the dark matter power spectrum as the figure of merit. The first 
method 113 811 is based on the theoretical prediction of the rela- 
tion between the halo mass and the bias factor (e.g., lf70l I88T - 
l9lll ). While the large-scale clustering measurements provide 
only the relative bias, galaxy-galaxy lensing measurements on 
small scales can provide robust estimates of the mean mass of 
the galaxy samples. The theoretical prediction of galaxy bias 
can in principle be refined to arbitrary accuracy by using nu- 
merical simulations, and its predictions depend on cosmology. 
Comparison of its prediction based on the lensing measure- 
ments of mass to the clustering measurements then determines 
the acceptable sets of cosmological parameters. This method 
works best at relatively high masses corresponding to the LRG 
sample or small clusters, at which the mass-dependence of 
the bias function on cosmology is strong and the uncertainties 
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in mass measurements are smallest due to their large volume 
coverage and high mass. Moreover, this method can be easily 
extended to the future surveys without much difficulty in the- 
ory and observation. The main observational systematic issue 
related to the method are the centroiding issue, which can be 
ameliorated if one uses lensing statistics that do not use small 
scale lensing information 1 20] . On the other hand, we have 
shown that this method suffers less from the overall lensing 
calibration error than the other weak lensing based methods. 
The method also assumes that the halo bias depends on the 
halo mass only. This is true by definition if all the halos in a 
certain mass range are included in the analysis, otherwise one 
must resort to halo selections based on an observable that is 
not sensitive to halo bias variations 17511 . 

The second method is to combine the large-scale measure- 
ments of galaxy-galaxy lensing and galaxy clusterin g, w hich 
yields the linear theory constraints on <7gf^ 55_1 125. 27j p76ll . 
In the formulation used here the method has very few theo- 
retical assumptions, since it is based on large scale clustering 
and lensing, where the two probes trace the same LSS with 
the cross-correlation coefficient close to unity [27]. While the 
measurement uncertainties depend on mass via its volume and 
total abundance, this method works well for a broad range 
of mass, with the best constraints in SDSS corresponding to 
the range between the L* galaxies to the low mass clusters, 
which allows multiple galaxy samples to be combined to sta- 
tistically tighten the resulting constraints. Furthermore, this 
method can be greatly improved with larger number density 
of background source galaxies at higher redshift in the future 
surveys, as its measurement uncertainties lag behind the clus- 
tering measurement uncertainties in the SDSS. In contrast, 
the cosmological constraint is the combination of <J%VL m , such 
that other constraints should be combined to break the degen- 
eracy. The future galaxy surveys we considered improve the 
constraints derived for the SDSS, mainly by going deeper in 
redshift and covering larger sky. However, the photometric 
galaxy surveys like the DES lack capability to map galaxies 
in three-dimensional positions and measure the angular clus- 
tering measurements, which become the limiting factor given 
the substantial improvements in lensing measurements. 

The third method is to utilize the (galaxy or cluster) abun- 
dance information, in conjunction with the lensing measure- 
ments of mean masses of galaxy samples. Due to the expo- 
nential sensitivity of the halo mass function to cosmological 
parameters in the clusters range, this method can yield some 
of the best constraints among the three methods, and the mass 
range of the LRGs to low mass clusters works best for this 
method too. The constraints are, however, comparable to the 
combined method I+II, which use both large-scale and small- 
scale lensing information, in addition to the large-scale clus- 
tering information. While the abundance is modeled by as- 
suming a log-normal scatter in the mass-observable relation, 



our investigation of the systematic errors in this method shows 
that the systematic errors such as the presence a g of invisi- 
ble halos and the skewness S3 of the probability distribution 
in the mass-observable relation are relatively innocuous given 
the present level of uncertainties in those parameters. For the 
future surveys this method continues to yield competitive con- 
straints, but is not superior to the combined method I+II. The 
methods I and II do not use abundance information and are 
thus not sensitive to the systematics associated with that, such 
as identifying multiple clusters in the same halo. 

We also compared the derived constraints to those from 
cosmic shear measurements that use shear-shear correlations. 
These can use information to significantly smaller scales and 
as a result can provide tighter constraints, despite the two- 
dimensional analysis assumed here. However, we argue that 
some of the gains cannot be achieved because of the bary- 
onic effects [84], and more importantly, the sample variance 
and the nonlinear evolution of the matter trispectrum make 
the lensing power spectra on small scales highly correlated 
rt8"oTl82"tl . setting the lower limit to the figure of merit one can 
obtain from the lensing power spectrum. Consequently, the 
constraints from the cosmic shear measurements are as strong 
as that from the combined method I+II. Furthermore, shear- 
shear analysis may also be more sensitive to the various spuri- 
ous systematics such as variable PSF Considering the serious 
systematic uncertainties present in the cosmic shear method it 
is useful to consider alternative methods, and the three meth- 
ods explored in this paper are the best alternatives proposed 
so far in measuring perturbations of dark matter, apart from 
redshift-space distortions which we do not consider in this pa- 
per. Their combined power is comparable to that of the shear- 
shear power spectrum, with very different systematics errors. 
In particular, the combined lensing and clustering analysis 
method, which combines the methods I and II of this paper, 
is a method that has not been discussed much in the literature, 
yet it gives predicted errors for SDSS that are comparable to 
the better known cluster abundance method with weak lens- 
ing calibration (method III). The method is currently being 
applied to the SDSS data and we expect this method will play 
an important role in the future surveys as well. 
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